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Abstract 

Microscopical approach to a deuteron breakup at high and intermediate ener- 
gies is proposed. We show that the quark exchange effects, resulting from the full 
asymmetry of the 6g-deuteron wave function with respect to the pair permuta- 
tions of quark variables, strongly affect the proton momentum distribution in the 
deuteron, as well as the polarization observables of inclusive deuteron breakup, 
when the "internal momentum" in the deuteron is of order of a few hundreds 
MeV/c. 



1 Introduction 



During the last decade the breakup of relativistic nonpolarized and polarized deuterons 
on hydrogen/nuclear target have been a subject of experimental and theoretical studies. 
They, as well as studies of e<i-scattering at high Q 2 , bring light on many aspects of the 
nucleon-nucleon interaction and the deuteron structure at short distances, including 
relativistic effects, effects of non-nucleon degrees of freedom, etc. The deuteron breakup 
experiments, we are discussing in this paper, include measurements of differential cross- 

section, E p ^0-, and polarization observables (the tensor analyzing power, T 20 , and the 
polarization transfer coefficient, Kq) of the (d, p) reaction with the final proton registered 
at zero angle and the momentum p > x (deuteron lab. momentum) at deuteron 
beam energy of a few GeV. 

We think that there are at least three kinematical regions where different physical 
pictures determine the reaction. At the region of the "internal momentum" in the 
deuteron, k < 0.1 GeV/c, the relevant degrees of freedom in the deuteron are nucleons 
and mesons. Here the nonrelativistic impulse approximation (IA) improved by the 
rescattering effects and the final state interactions, both calculated with the standard 
two-nucleon deuteron wave functions, seems to be valid. At the region of very high k 
(according to estimations of k > 1 GeV/c) one has to provide a smooth transition 
to the perturbative QCD. 

In the present paper we discuss the deuteron breakup at the intermediate region 
(k between 0.1 and 1 GeV/c) where the nucleons in the deuteron are assumed to loose 
gradually their individuality, but the pure perturbative QCD does not yet work. At 
this region the deuteron structure can be described in terms of the constituent quarks 
and the chiral mesons within the dynamical scheme similar to that used in the baryon 
spectroscopy [0. Then one should note that the Pauli principle at the constituent 
quark level gives rise to a number of the short-range baryon-baryon components (AW*, 
N*N, N*N*) of the deuteron wave function [3, HI. 

However, instead of the usual Fock probabilities for the baryon-baryon components, 
which are applicable only for the structureless (or well separated baryons), one should 
use a language of "effective numbers" [|J developed in nuclear cluster physics ||. 

In this work we shall restrict ourselves to a simplified model, in which the quark 
antisymmetrization effects are taken into account for the deuteron S'-wave only. Due 
to the centrifugal effect the quark-exchange contribution for the .D-wave is estimated 
to be three orders of magnitude less than the one for the .S'-wave ||. Besides, we do 
not consider here a delicate question of the relativistic internal momentum k in the 
deuteron: we use it in the same way as for the equal-mass NN configuration in the 
framework of the light-cone dynamics ("minimal relativistic prescription"). The 3gx3g 
decomposition of the deuteron six-quark wave function contains, apart from the equal- 
mass NN configuration, a set of configurations, NN*, whose "constituents" (N and 
N*, respectively) have different masses. Meanwhile, the main result of this work is that 
the quark exchanges between 3q clusters in the deuteron give a visible contribution in 
the differential cross-section and polarization observables of the (d,p) breakup in the 
intermediate region, and have to be taken into consideration. 



The outline of the paper is organized as follows. In Sec. 2 we consider a decompo- 
sition of the deuteron six-quark wave function into 3qx3q clusters and how it modifies 
a matrix element of the deuteron breakup. In Sec. 3 expressions for the proton mo- 
mentum distribution in the deuteron and polarization observables are obtained. The 
results of numerical calculations and conclusions are given in Sec. 4. 




Fig. 1. Deuteron inclusive breakup in spectator approximation. 

2 NN and NN* components of the deuteron six- 
quark wave function 

The deuteron wave function, considered at small internucleon distances as a six-quark 
(6g) object, was shown || |7|, |] to be qualitatively equivalent to the Resonating Group 
Method (RGM) wave function 

2, . . . , 6) = A {<p N (l, 2, 3)^(4, 5, 6) X (r)} , (2.1) 

where A = — 9P 36 ) is a quark antisymmetrizer and ^(1,2,3) and y>jv(4, 5,6) 

are the wave functions of the nucleon three quark (3g) clusters; x( r ) is the RGM 
distribution function. 

The presence of the quark antisymmetrizer A in expression ( |2.1| ) is the main dif- 
ference between the microscopical quark and the meson-nucleon points of view on the 
deuteron structure at short distances. In the absence of the quark antisymmetrizer the 
RGM distribution function would coincide, up to the renormalization effects at short 
range, with the conventional deuteron wave function. 

Due to the antisymmetrizer the deuteron wave function ( |2.1|) , being decomposed 
into 3q x 3q clusters, includes, apart from the standard pn component, the nontrivial 
NN*, N*N and N*N* components which respond to all possible nucleon resonance 
states, N* (see, e.g., [§J). So, in the spectator approximation (Fig.l) the reaction 
matrix element is 

M dp ^ p r X = s/2 £ i> d Bp ,(k p ,)r Bp ^ x , (2.2) 

B,M B 

where the summation over B includes pn and all pN* configurations of the deuteron; 
Mb is ^-projection of 5-baryon spin. The other notations in Q2.2p are as follows: 



r_B P ^x is the amplitude of the reaction B + p — > X: 

r P ^x = / n n rf3 ^ rf3 ^- e-^^+v-^)^ 

i=lj=7 

x 2, 3, 7, 8, 9)|f(l, 2, 3, 7, 8, 9)|p B (l, 2 , 3)p p (7, 8, 9)), (2.3) 

rx = igr, + gr.). (2.4, 

and <^b(1, 2, 3) is a wave function of the 3q cluster of the baryon B; <^tv(7, 8, 9) is the 
same for the target proton; stands for the i-th quark coordinate and 

V4>'(V) = 7^372 / d " r e ' v >V( r )> ( 2 - 5 ) 

is a Fourier transformation of the overlap, ^^(r), between the 6q wave function ( |2.1| ) 
of the deuteron and the pB wave function; r stands for the relative coordinate between 
the target and the baryon B: 



6! 



1/2 



</(r)= I^J (^(1,2,3)^(4,5,6)|^(1,2,3,4,5,6)), (2.6) 
r = X>-X>, (2.7) 

i=l j=4 

In [|J the expressions for the overlaps ^^ p /(k p /) were obtained for all pB configurations 
produced by the antisymmetrization of the S'-wave component of the deuteron wave 
function. As the first step to an appropriate choice of the RGM distribution function, 
x(r), we use the conventional NN deuteron wave function, Xnn(y), consisting of the 
S- and D- components, and modify it according to the standard RGM renormalization 
condition || |T0 



A^r.rOWrVV, (2i 



where N(r, r') - is a norm operator 

A>(r,r') = 5(r- 0-9(^(1, 2, 3)^(4, 5, 6)} 5 T=io5(r-r")| 

xP 3G \{<p N {l, 2, 3)^(4, 5, 6)} ST=10 5(r' - r")). (2.9) 

Fig. 2 displays the ratios X L=0 ( r ) / Xnn ( r ) an d X L=2 ( r ) / Xnn ( r )- I* 1 the further cal- 
culations we approximate the radial parts of the S- and .D-components of the RGM 
distribution function by gaussians: 

^^ = 5> fc e- Q < ^^ = r 2 ]T5 fc e-^ 2 . (2.10) 
r k r k 

As the nucleon and excited states wave functions we have used the wave functions of 
the Translationally Invariant Shell Model JTTJ for three particles, i.e. the harmonic os- 
cillator wave functions with oscillator parameter b, with the center-of-mass oscillations 



being removed. All detailes can be found in M . The parameters of the approximation 
( |2.10 ) for the RGM distribution function based on the Paris potential are given in 



Table 1. Making use the gaussian approximation all necessary calculations could be 
performed analytically jp and the Fourier transformation Q2.5Q is given by 

V4 P '(kp') = x(K')$b,n + 3 E (L B S B m B fi B \ J B M B ) 



m B ^B ^dim[f B ] 
1b E A k I^ BLB , 00 (k p r,a k ) n sms (V) (-i)i+s B+ zr B 



x^f(2T B + 1)(2S B + 1){S b -VbV p > I lM d ) 



x E {[fB]S B TB\[2}S l2 T l2 - l -\) 

Sl2=Tl2=0,l 



545=T 45 =0,1 

x( _i)S 12 +s 45 , > ■ „ 2 , . jS - , i , (2.u: 




For the explanation of notations in ( 2.11 ), see Q. The functions I^ B i o are sum 

B x B ' 



marized in Appendix A. The fractional parentage coefficients (CFP), 7^, for N B = 3 
could be found in Appendix B (the CFP for N B = 0, 1, 2 are given in @]).The physical 
meaning of the parameter b appearing in these functions is the m.s. radius of a nucleon 
quark core, i.e. a characteristic distance where quark exchanges becomes noticeable. 
In Table 2 we summarize the baryons B which are taken into account in our calcula- 
tions and give the effective numbers, N^ p , for three values of the parameter b; for the 
definition of the effective number N^ p see, e.g. ||. The baryons incorporated in the 
components with the effective numbers N^ p < 10~ 5 were ignored. 

It should be mentioned here that the baryons from lines 2-5 and 7-10 of Table 2, 
have the negative parity and thus produce effective P-waves of the relative motion in 
the deuteron. It was already mentioned (|L2| and [|13|]) that the P-waves should strongly 



affect polarization phenomena of the reaction at the region of internal momentum in 
the deuteron of few hundreds MeV/c. 

There are at least two problems connected to the relativistic description which we 
would like to discuss here shortly. First, RGM is based on nonrelativistic dynamical 
equations. To establish a correspondence between the observed momentum of the spec- 
tator and the wave function argument we use the "minimal relativization prescription" . 
This means that we identify the light-cone variable k with the relative momentum of 
the 3q clusters in the deuteron, considering their dynamics nonrelativistically, but in 
terms of new variable k (see, e.g., 01). Second, the different components of the 
deuteron wave function (NN, NN* etc.) consists of "particles" (3q clusters) with dif- 
ferent masses. In this case every component, involved in our calculations, in the limit 
of infinite distances should depend on its own "internal momentum". However such a 
limit is allowed only for the usual n-p component and the effective NN* components 
appear only in the nucleon overlap region and they are direct consequence of the under- 
lying quark structure of a baryon. Thus this question cannot be solved at the baryon 



level. In the present work we assume for simplicity that the " internal momenta" in the 
n-p and iV — iV* chanels are the same. 



3 Momentum distribution of the protons in the 
deuteron and polarization observables in the 
spectator approximation 

According to the spectator approximation Q2.2| ) the reaction matrix element squared is 
J2 \M dp ^ p , x \ 2 = 2a J2 |V4/(V)f • ( 3 - 12 ) 

X B 

In ( [3.12j ) it is assumed that 

^ FpB'-^x FpB^x = &o5bb>, (3.13) 
x 

where o"o is of order of the total NN cross section. This means that the size of B is 
supposed to be the same as the proton's one and inelasticities in baryon-baryon colli- 
sions are small at the energy scale of few GeV. Now the proton momentum distribution 
in the deuteron is given by 

n{k p/ ) = u 2 + w 2 + u\ + v\ + ~v\ + v\ + ~v\ + 1v\ + v\ + 2v 2 + vj, (3.14) 

where 

u = X L= \K>) A k / °o,o (V; a k ), (3.15) 

w = x L=2 (£y), (3-16) 

4 

«1 = jT/= E A k J 20,Oo( fc p'; a fc)' ( 3 - 17 ) 

4 

vi = v 2 = v$ = v 6 = -J2 A k Ai,oo(V; a k)i (3-18) 
y k 

2 

u 3 = v 4 = v 7 = v s = -VsJ2 A k ^31,00 (V; a *)- ( 3 - 19 ) 

y fc 

The differential cross section of the reaction is proportional to the proton momentum 
distribution, n(k p i), (see, e.g., and the tensor analyzing power, T 2 q, and polariza- 



tion transfer coefficient, kq, are given by 

rp 1 gx/2gg - w 2 - H - |yj + \v 2 + |^ 8 2 

J 20 — — 7= 71 s > (^.ZUJ 

K = (W 2 - U> 2 - ^MW + u\ + f 2 - 2 + W 2 - -U* + 2W 2 

- ^ 2 + 2t; 2 - ±-vl)/n(k pl ). (3.21) 



4 Numerical calculations and conclusions 



The results of our calculations and comparison with the experimental data are shown 
in Figs. 3-9. Figs. 5 and 7 explore the dependence of the polarization observables T 20 
and k,q on the parameter b of the oscillator quark potential. In Figs. 3, 4, 6 and 8 we 
demonstrate insensitivity of the current approach on the choice of input NN potential. 

In Fig. 9 we compare the results with the experimental data on the T 20 — k plane. 
It was recently shown [12], that a correlation between T 20 and Kq should give important 
information about the non-nucleon degrees of freedom in the deuteron. Particularly, 
in the case, when the deuteron wave function consists of the two standard (S— and 
D— wave) components only, the parametric curve describing the correlation between 
T 20 and k,q must lie on the circle. Our calculations give an example of a deformation 
of this relation produced by additional components of the deuteron wave function. 

The points for the proton momentum distribution in the deuteron n(k p >) (Figs. 3 
and 4) were extracted from data for the p(d,p) reaction cross section Jl5[ (see also this 
reference, as well as flSfl , fT7|| for A(d,p) data). T 20 was measured in [jP7||-p2"| and 



m 



21, 23, 24, 25 



The conclusions of the present work follows: 

(i) When the internal momentum in the deuteron is of a few hundreds MeV/c, the 
effects of the quark exchange between three-quark clusters in the deuteron are as im- 
portant as the ones originating from the relative motion of the nucleons (estimated by 
the standard two-nucleon potentials). 

(ii) The pN* components produce effective P-waves in the deuteron which correct the 
momentum distribution of protons in the deuteron and polarization observables of the 
deuteron breakup at high/ intermediate energy in the "right" direction. 

(iii) The results of the current approach depend weekly on the choice of an input two- 
nucleon potential, but, at the same time, very sensitive to the quark core radius of 
the nucleon, b. However, it should be noted the value b = 0.8 fm which provide the 
best fit to data is somewhat larger than the m.s.r. of the quark core of a free nucleon, 
b = 0.5 4- 0.6 fm, commonly used in theoretical calculations. 

It should be also stressed that we consider here the intermediate region between 
the pure NN and the perturbative QCD regimes and our model cannot be used for a 
region of very high k. 
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Appendix A 



/ 27 N§ f W 15 + W 

Voo " Us + Wa k w) beXP \ 12 3 + 5a k W 

^11,00 = -^00,00 9i b k 



1 20,00 ~~ 1 00,00 



n°nnn V6 ( g - 3 -g l9 + 4bW 



2 J J 2 



f 2 - T° 4 -n 2 b 2 k 2 
1 22,00 ~~ J 00,00 \/ 



4,00 = ^o°o,oo V^i b k (g — lg x g + M^fc 2 ) , 



where 



4a,6 2 - 3 4a,6 2 - 3 



9 = TT—Tr — T5' 9i = 



15 + 16a fc 6 2 ' y 18 + 30a fc 6 2 ' 



Appendix B. 

Baryon orbital wave functions |3(A//)[/]L(r)) 
The choice of Jacobi coordinates is 



ri - r 2 ri + r 2 - 2r 



r l2 = 7S— , P 



3 



V2 ' " v/6 



|3(30)[3]l(lll)) = ^& 2 oo(ri 2 )¥'iiAf(p) + i £ (21 | 1M> 

^ V o V3 mi , m2 

xy?22 mi (ri2)y?iim 2 (p) - ^ooo(ri 2 )v^3iM(p) 
|3(30)[21]l(112)) = -^2oo(ri2)^iiM(p)-i £ (21 m x m 2 \ IM) 

^ " mi,rri2 

xv?22 mi (ri2)v ? iim 2 (p) - ^-V ? ooo(ri 2 )v ? 3iM(p) 



|3(30)[21]1(121)> 



|3(11)[21]1(112)) 
|3(11)[21]1(121)> 
|3(30)[1 3 ]1(123)) 

|3(11)[21]2(112)> 
|3(11)[21]2(121)> 

|3(30)[3]3(111)) 



y/5 

-^-</?31M(ri 2 )</?OOo(p) frfUM (ri2)V200 (p) 

i (12 miffl 2 1 lM)</? llmi (r 12 )</? 2 2m 2 (p) 



mi ,m2 



2 v^5 

^<^20o(ri2)v ? iiAf(p) o~ XI ( 21 m l m 2 I 1M) 

" " mi,m2 

X^22mi(ri 2 )^llm 2 (p) 

2 

^V 5 iiM(ri 2 )v^20o(p) ~- X ( 12 mim 2 | 1M) 

X</?llmi(ri2)^22m 2 (p) 

1 l/5" 

-y?3iAf(ri2)y?ooo(p) - 2Y gV ? iiM(ri2)^20o(p) 
— 7= X (12 mim 2 I lM)v?ii mi (ri 2 )v?22 m2 (p) 



- J] (21 mps I 2M)v? 22mi (r 12 )v? llm2 (p) 

mi,m2 

X (12 mim 2 I 2M)v?ii mi (ri 2 )v? 22m2 (p) 



mi,m2 

V3 



X (21 mim 2 I 3M)v? 22mi (ri 2 )v?ii m2 (p) 



mi,m2 



--y?ooo(ri 2 )y?33Af(p) 



|3(30)[21]3(112)) = -- E (21m 1 m 2 \3M)ip 22mi {r 12 )cp llm2 (p) 



rV ? ooo(i"i 2 )^33M(p) 



|3(30)[21]3(121)) 



|3(30)[1 3 ]3(123)) = 



-y ^V 5 33m(i - 12)v ;, OOo(p) 

~ (12 m 1 m 2 I 3M)v9 llmi (r 12 )v? 22m2 (p) 
^33Af(ri 2 )v?ooo(p) 

— 2~ Yj (12 mim 2 | 3M)</? llmi (r 12 )</? 22m2 (p) 



Table 1. 



The parameters of the approximation (|2.10| ) for the RGM distribution 
function based on the Paris potential. The oscillator parameter b = 0.8 fm. 



n 




ot n {fm 2 ) 


B n (fm-y' 2 ) 


Pn (fm- 2 ) 


1 


1.3863 • 1Q-' 2 


8.845 • lO" 3 


3.1946 ■ 10~ 4 


4.1231 ■ lo" 2 


2 


7.1097- 10" 2 


3.1668- lo" 2 


3.6241 • 10" 3 


1.1980- 10- 1 


3 


1.7493 ■ lO" 1 


1.0619- 10- 1 


1.7823 ■ lo" 2 


2.6855- 10- 1 


4 


2.9944 • 10" 1 


3.8459- 10- 1 


7.5863 • lo" 2 


5.5486- 10- 1 


5 


-6.8075 • 10" 1 


2.4318 


3.3025 • 10" 1 


1.3350 


6 


1.8345 • 10- 1 


3.5360 


-2.6708 • 10" 1 


5.4722 



Table 2. 

The effective numbers N^ p of different baryon-proton configurations in the 
deuteron. The calculations were done for the Paris potential. 





J P T 


State 




b = 0.5 /m 


b = 0.7 fm 


b = 0.8 fm 


1 


1+ 1 
2 2 


|0(00)[3]oii) 


0.982 


0.996 


1.005 


2 


3- 1 
2 2 


H(10)P1]1~> 


1.52 x 10- 3 


3.60 x 10- 3 


4.46 x 10- 3 


3 


1- 1 
2 2 


|i(io)[2i]iii) 


7.59 x 10~ 4 


1.80 x 10- 3 


2.23 x 10- 3 


4 


1- 1 
2 2 


, / NT , 31. 

H(10)[21]l--) 


3.79 x 10" 4 


9.01 x 10" 4 


1.16 x 10- 3 


5 


3- 1 
2 2 


, / NT , 31. 

H(10)[21]l--) 


7.59 x 10~ 4 


1.80 x 10- 3 


2.23 x 10- 3 


6 


1+ 1 
2 2 


|2(20)[21]o||> 


2.32 x 10- 3 


5.47 x 10- 3 


6.75 x 10- 3 


7 


1- 1 
2 2 


|3(30)[21]l||> 


2.52 x 10" 4 


3.63 x 10~ 4 


3.67 x 10" 4 


8 


3- 1 
2 2 


|3(30)[21]1~> 


5.04 x 10~ 4 


7.25 x 10~ 4 


7.33 x 10- 4 


9 


1- 1 
2 2 


, / u 1 31, 
|3(30)[21]l--> 


1.26 x 10" 4 


1.81 x 10 -4 


1.83 x 10" 4 


10 


3- 1 
2 2 


, / u 1 31, 
|3(30)[21]l--> 


2.52 x 10" 4 


3.63 x 10" 4 


3.67 x 1Q- 4 
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Fig. 2. Ratios \ L C \ r ) l 'Xn n '( r ) and x L 2 ( r )/XAw ( r )- Calculations are performed with 
the Paris potential and the parameter b = 0.8 fm. 



Fig. 3. Momentum distribution of the proton in the deuteron. Solid line stands for 
calculations within the current approach with b = 0.8 fm, dashed and short-dashed 
lines - to b = 0.5 and 0.7 fm, respectively. Calculations are performed with the Paris 
potential. The points are extracted from l H{d 1 p) cross section data [O . 



Fig. 4. Momentum distribution of the proton in the deuteron. Solid and dashed 
lines stand for calculations within the current approach with b = 0.8 fm, with Paris 
and Reid- soft- core potentials, respectively. By short-dashed line we denote the curve 
obtained in the IA with Paris potential. 



Fig. 5. T 2 o in the current approach. Solid line corresponds to b = 0.8 fm, dashed and 
short-dashed lines - to b = 0.7 and 0.5 fm, respectively. Calculations are performed 
with the Paris potential. The experimental data are from [|T7|]-||T9 • 



Fig. 6. T20 in the current approach (b = 0.8 fm), with the Paris potential (solid line) 
and the Reid-soft-core potential (short-dashed line). Dashed line stands for the IA 



calculations with the Paris potential. The experimental data are from [|T7|-||T9" 



Fig. 7. The same as in Fig. 5, for kq. The experimental data are from pi, 53, 21 



Fig. 8. The same as in Fig. 6, for kq. 



Fig. 9. The k,q vs T 20 plot in the current approach (dashed line). The solid line 
represents the results of IA ||12|| . Calculations are performed with the Paris potential. 
Data are from ||23|| . 
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